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VARIATIONAL CALCULUS FOR HYPERSURFACE FUNCTIONALS: 
SINGULAR YAMABE PROBLEM WILLMORE ENERGIES 

MICHAEL GLAROS**, A. ROD COVER* 1 , MATTHEW HALBASCH** & ANDREW WALDRON** 


Abstract. We develop the calculus for hypersurface variations based on variation of 
the hypersurface defining function. This is used to show that the functional gradient 
of a new Willmore-like, conformal hypersurface energy agrees exactly with the ob¬ 
struction to smoothly solving the singular Yamabe problem for conformally compact 
four-manifolds. We give explicit hypersurface formulae for both the energy functional 
and the obstruction. 


Keywords: Calculus of variations, conformally compact, conformal geometry, hypersurfaces, Willmore energy, 
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1. Introduction 

Hypersurface geometry is critically important for the analysis of boundary problems 
in mathematics and physics. Conformal hypersurfaces form a subtle but important class 
of cases, particularly because of their role in treating the boundary at infinity of confor¬ 
mally compact structures. Naturally the fundamental problem is the construction and 
study of local and global invariants. The work |GW15j develops an effective approach to 
hypersurface geometries by treating them as conformal infinities (see also the announce¬ 
ment JGW131 ). This approach includes a calculus for hypersurface invariants based on 
defining functions. 

In this article we develop the corresponding variational theory; namely a hypersurface 
variational calculus based around defining functions. A main result is Theorem 13. 11 1 which 
provides the general formula for computing variations of embeddings of hypersurface 
functionals by varying the underlying hypersurface defining function. A key motivation 
for our development of that theory was to develop an effective approach to computing the 
Euler-Lagrange equations corresponding to the new conformal energies of |GW15j . These 
functionals generalise to higher dimensions the Willmore energy |Wil65| . also known as 
the rigid string action [Pol86 ]. For the crucial case of four manifolds, we show that the 
gradient of this functional agrees with the obstruction to smoothly solving the Yamabe 
problem for conformally compact structures. 

The Willmore energy and its associated gradient provide important invariants for sur¬ 
faces. Indeed the Willmore conjecture concerning absolute minimizers of this energy has 
stimulated much geometric analysis |Riv08| : the general proof of the conjecture was pro¬ 
vided by Marques-Neves |MN14l fMN13| . In physics, the Willmore energy functional has 
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applications in string theory where it appears as the entanglement entropy for a three 
dimensional conformal field theory |AGS14| . It also arises as a certain log coefficient 
in a holographic notion of physical observables, and is closely related to the so-called 
conformal anomaly |GW99j . Furthermore it appears as the renormalized area of a min¬ 
imal surface embedded in a hyperbolic three-manifold |AM10| . An important aspect of 
these invariants is their conformal invariance. While the Willmore energy integrand is 
not surprising its gradient, with respect to variations of embeddings, is an important 
local conformal invariant because it has linear leading term. We call this gradient the 
Willmore invariant. 

These applications and observations suggest a potentially significant role in geom¬ 
etry and physics for corresponding conformal invariants in higher dimensions. For four¬ 
dimensional hypersurfaces in Euclidean spaces, conformal energy functionals have been 
provided in |Guv05| and |Vyal3|. Recently it was shown that for hypersurfaces of every 
dimension, there exists a conformally invariant canonical analog of both the Willmore 
energy and the Willmore invariant. These are generated and constructed uniformly via 
a singular Yamabe problem that places the hypersurface as the boundary at infinity of a 
conformally compact, constant scalar curvature, manifold (GW1311GW15) . In each higher 
dimension, the analog of the Willmore invariant is called the obstruction density. This 
work was partly inspired by the results [ACF92] of Andersson-Chrusciel-Friedrich who 
studied boundary asymptotics for this Loewner-Nirenberg-type problem |LN74j . They 
found that the three-manifold obstruction to smoothness was a conformal invariant and 
also gave some information about the corresponding obstructions in higher dimensions. 
For three-manifold boundaries the obstruction density is easily shown to agree with the 
Willmore invariant. 

In |GW131 IGW15] the singular Yamabe problem was recast as the tool for treating 
hypersurface conformal geometry that is analogous to the Fefferman-Graham Poincare- 
metric as an approach to intrinsic conformal geometry. From this perspective, and using 
tractor calculus to re-express the Yamabe equation, the problem of computing the ob¬ 
struction is reduced to a simple algorithm that also reveals the qualitative aspects of the 
obstruction densities and higher Willmore energies. In particular there is a dimension 
parity dichotomy: For each even hypersurface dimension, the obstruction density is a 
linear-leading-order conformal invariant that is a close analog of the Willmore invariant 
(and should be considered as a fundamental scalar curvature invariant); via elementary 
representation theory, odd dimensional hypersurfaces cannot admit linear leading order 
invariants. 

It is an interesting problem to understand the higher Willmore energies and obstruction 
densities for odd dimensional hypersurfaces. A particularly important case for both 
mathematics and physics is when these are four-manifold boundaries. In this article we 
employ the general holographic formula for the obstruction density | IGW15| to compute 
the four-manifold obstruction density: 


Proposition 1.1. The obstruction density for d = 3 is given by 


B, = - 


6 L 


0 ™ o , fl L o 0 7T- n ° 7 ° J 7T. fi 

L c ^Kb) 0 -W ab ) ~^ ab B ab +K 2 -7W ab nl b +2W ab W ab +R ab W d W cabd +Wj bc W' 


abc 


Theorem 15.101 recalls that is a conformally invariant density of weight — d, and our 
tensor notations are explained in Section 11.11 below, while the hypersurface invariants 
appearing in the above formula are introduced in Section [2j In particular II denotes the 
trace-free second fundamental form, while L ab is a second order, conformally invariant 
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differential operator (that appears in Bernstein-Gerfand-Gel’fand (BGG) sequences) 
such that L al) ll a b = V“V b II a 6 + lower order terms. 

As mentioned above, another main result in |GW15| is the construction of conformally 
invariant functionals that generalise the Willmore energy to higher dimensions. For a 
hypersurface £ of dimension d these take the form 

2 3 = a a [ dA 3 N A P j N A , 

J s 

where aj is a non-vanishing constant and Pj is a uniformly constructed family of con¬ 
formally invariant operators, and N A is the normal tractor of [ BEG94 J introduced in 

Section [5.21 For d even, the operators Pj take the form (A T )z + lower order terms , so 
are hypersurface analogs of the GJMS operators of Graham, Jennes, Mason and Spar¬ 
ling |GJMS92) . For d both even and odd a holographic formula for Pj is known. This 
gives that P 3 takes the form Il ab VjV^ + lower order terms (the detailed formula appears 
in Equation (16.21) 1. So for generic hypersurfaces P 3 is Laplacian-like. Using these results 
for P 3 yields 

( 1 . 1 ) = dA sA ah T ah , 

where the conformally invariant Fialkow tensor takes the form T a \, = IP. + lower 
order terms ; see Equation (12. 7[) and 03 = | has been chosen for later convenience. 

Our final main result is that the functional gradient of the energy Z 3 agrees precisely 
with the obstruction £>3 to the singular Yamabe problem mentioned above: 

Proposition 1.2. With respect to hypersurface variations, the gradient 0 /Z 3 is B 3 . 

Our article is structured as follows: In Section [2j we review basic hypersurface theory 
and use conformal invariance as the organising principle when producing expressions for 
hypersurface invariants (including a novel hypersurface Bach tensor B a f, which already 
appears in Proposition 11.11 above!. In Section [ 3 ] we develop our hypersurface variational 
calculus and treat Willmore and minimal surfaces as basic examples. In Section [I] we 
compute the variational gradient for our four-nranifold analog of the Willmore energy. 
Section [5] is devoted to a review of the essential tractor calculus methods needed to ef¬ 
fectively treat the singular Yamabe problem. In Section [ 6 ] we compute the obstruction 
density for four-manifold hypersurfaces. We also compute the obstruction for example 
metrics, including a hypersurface surrounding the horizon of a Schwarzchild black hole. 
The appendix gives explicit formulae for the obstruction density for manifolds of dimen¬ 
sions five or less (so including £> 4 ) that can easily be implemented for explicit metrics 
using computer software packages. 

1.1. Riemannian and conformal geometry conventions. We work with manifolds 
M of dimension d and embedded hypersurfaces of dimension d := d — 1. When the 
dimension d equals three or four, we often refer to the latter as surfaces and spaces, 
respectively. (Note that the exterior derivative will be denoted by d, to avoid confusion 
with the dimension d.) When M is equipped with a Riemannian metric g, its Levi-Civita 
connection will be denoted by V or V a in an abstract index notation (c/. [PR84| ). The 
corresponding Riemann curvature tensor R is 

R(u, v)w = [V u , V v \w - V [ U>V ]W , 

for arbitrary vector fields u, v and w. In the abstract index notation, R is denoted by 
R a b c d and R(u,v)w is u a v h R a i> c dW d . Cotangent and tangent spaces will be canonically 
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identified using the metric tensor, g a b and this will be used to raise and lower indices in 
the standard fashion. 

The Riemann curvature can be decomposed into the trace-free Weyl curvature W a b c d 
and the symmetric Schouten tensor P a b according to 

Rabcd Wabcd T 2(/ a[cRd]b ^9b[cR d]a ■ 

Here antisymmetrization over a pair of indices is denoted by square brackets so that 
X[ a b] ■= ^(X a b — Xba). The Schouten and Ricci tensors are related by 

Ric bd ■= Rab a d = (d - 2)P bd + Sab J , J := Pa • 

The scalar curvature Sc = g ab R\c a b, thus J = Sc/(2 {d — 1)). In two dimensions the 
Schouten tensor is ill-defined, but we take J = ^ Sc in that case. In dimension d > 3, the 
curl of the Schouten tensor gives the Cotton tensor 

Cabc := 2V[(jPfrJc • 

In dimensions d > 4, 

(1.2) (d-3)C abc = V d Wab d c- 

A conformal geometry is a manifold equipped with a conformal structure c, namely 
an equivalence class of metrics with equivalence relation g' ~ g when g' = £l 2 g for 

C°°M 3 > 0. The corresponding Weyl tensors obey W° a b c d = I Vab c d and thus the 

Weyl tensor is a conformal invariant. Conformal invariants may also be density-valued A 

w 

weight w density is a section of the oriented line bundle [(A“TM) 2 ] 2d =: £M[w\. Vector 
bundles VM tensored with this are denoted VM <S> £M[w\ =: VM[w\. Each metric 
j£c determines a section of (f\ d T*M ) 2 and so there is a tautological global section g 
of Q 2 T*M[ 2] called the conformal metric. Hence, each g € c is in 1 : 1 correspondence 
with a strictly positive section r £ £M [1] via g = t 2 g. Such sections r will be referred 
to as a true scale, and we will often present conformally invariant formulae in terms of 
a Rienrannian metric g by making an (arbitrary) choice of scale r. It will be clear from 
context where this has been done. 

In situations where we are dealing with tensors built from higher rank tensors by 
contractions with vectors or covectors into the first slot such as X(u,-) (or equivalently 

u 

u a X a b), we will often use the notation X a - In higher rank cases where more indices 

are missing, an alternating leftmost-rightmost contraction procedure has been followed, 

u 

so if X is rank 3, then Xb denotes X(u,-,u) or equivalently u a X a b c u c . Contractions 
across multiple indices of tensors with one another will be denoted by an obvious bracket 
notation, for example, if Y is rank two and X is rank three, X(Y, •) denotes the covector 
XabcY ab . Also, for rank two symmetric tensors S,T, we will have occasion to employ 
matrix notations such as S 2 b := S ac S£ and trST := S a bT ba . Unit weight symmetrization 
over groups of indices is denoted by round brackets, which we adorn with the symbol o 
when projecting onto the trace-free part of this, for example Aq a ^ 0 := ^(X a b + Xb a ) — 
2g a bX c c ', for mixed symmetry tensors, projection onto the trace-free part of a group of 
indices will be indicated by the notation { }o and their corresponding section spaces will 
be decorated with a subscript o. Finally, we will also use a dot product notation for 
inner products of vectors g(u,v) = u a v b = u.v as well as divergences V .Xb := X a X a b, 
and |«| := \Ju a u a := x/ri 2 denotes the length of a vector u. 
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2. Embedded hypersurfaces 

A hypersurface S is a smoothly embedded codimension 1 submanifold of a smooth 
manifold M. A function s € C°°M is called a defining function for £ if the hypersurface 
is given by its zero locus, i.e. £ = Z(s) and the exterior derivative ds =: n is nowhere 
vanishing along £. The covector held n|s is then a conorm.al to the hypersurface and 
h a := (u a /l n l)|v ^ ie unit conormal while n a \s is a norm.al vector. 

The conormal allows the tangent bundle T£ and the subbundle of TM |s orthogonal 
to n a , denoted TM , to be canonically identified. Thus we may use the same abstract 
indices for T£ as for TM. The projection of tensors, to hypersurface tensors will be 
denoted by a superscript T. Given an extension h a of the unit normal to M, we will also 
write v T := v — hh.v for vectors (and similarly for tensors) v € TM. Upon restriction, 
these give hypersurface vectors (and tensors). 

The metric g on M induces a metric g on £ given by 

(2-1) 9ab = 9abh ~ n a h b . 

Given a hypersurface vector held v a € T(T£), the Gauls formula relates the restriction 
to T£ of the Levi-Civita connection V on M to the hypersurface Levi-Civita connection 
of g according to 

(2.2) V a v b = Vjv b + n b H ac v c . 

In the above formula, the right hand side is defined independently of how v a € T£ is 
extended to v a € TM. In general, we will use the term tangential for operators O along £ 
if they have the property that for v a smooth extension of some tensor v dehned along £, 
the quantity Ov |s is independent of the choice of this extension. In the above display, 
the second fundamental form II a 5 € T(© 2 T*£) is given by 

Hab — n b ■ 

Its “averaged” trace is the mean curvature 


We will often use bars to distinguish intrinsic hypersurface £ quantities from their host 
space counterparts. In particular, the equations of Gaufi relate ambient curvatures to 
intrinsic ones according to 

^abcd hf,brd 2II r ," r , 

(2.3) Rl c a& = Hl c afc + Hafc — dHH ab , 

Sc — 2 Ric(n, h) = Sc — tr II 2 + d 2 H 2 . 

The last of these recovers Gauls’ Theorema Egregium. for a Euclidean ambient space. The 
Codazzi-Mainardi equation gives the covariant curl of the second fundamental form in 
terms of ambient curvature 

- 1 fl t 

(2.4) V [a n b]c = -R T cba . 

Finally the divergence of the Codazzi-Mainardi equation determines the Laplacian of 
mean curvature: 

- 1 - - _ h 

AH = -(V.V.n- V a RicI) . 

d 


(2.5) 
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2.1. Conformal hypersurfaces. The ambient conformal structure c induces a confor¬ 
mal structure c on X. Locally we may always assume there is a section h a € F(T*M[ 1]) 
obeying g ab n a n b = 1. This conformal unit conormal plays the role of a Riemannian unit 
conormal. The trace-free part of the second fundamental form Jl ab = II a fe — Hg abl defined 
(in a choice of scale) in terms of h a via Equation (12.21) , is a conformally invariant section 
of Q 2 T*M[ 1], Thus the quantity 

(2.6) K :=trfl 2 € T(££[—2]), 

is a conformal density of weight —2, that we call the rigidity density. 

Conformal invariance of the trace-free second fundamental form II and the Weyl tensor 
can be used to decompose the equations of Gau£ and Codazzi-Mainardi equations into 
conformally invariant parts. First, the Ricci equation in (12.31) gives what we shall call 
the Fialkow-Gaufi equation (c/. |Vyal3j) 

(2.7) n 2 ab - - Qab^zr ~ Wab = (d - 2 )( P ^ - Pab + HJ1 * b + 2 ~ gabH ) =: ~ 2 ^ ab ■ 

The left hand side of this equation is conformally invariant, which proves that the Fialkow 
tensor J- ab defined by the above in dimension d > 4 is a conformally invariant section 
of Q 2 T*M[ 0]. This combined with the trace-free second fundamental form gives a new 
distinguished density 

(2.8) L := tr(ILT) € T(£E[—3]), 

which we call the m.embrane rigidity density. 

In dimensions d > 4, the trace-free part of the first Gauft equation allows the tangential 
part of the ambient Weyl tensor to be traded for its hypersurface counterpart: 

I — o 2 / ® ^ ° 2 \ 2 r 

W a bccL = W a bcd ~ 2n a [ c II c (]; ) — - _ ^ (5a[eP-c(]& — g b [ c IL c ^ a ) + - _ _ 2 ) 9a[c9d]bK ■ 

In dimensions d > 4, the third Gaufi equation expresses the rigidity density in terms of 
the difference between ambient and hypersurface scalar curvatures: 

K = 2 (d - 1) (j - P(ft h) — J + ^ H 2 ^J . 

This shows that J — P(ft n) — J + ^ H 2 is a conformally invariant weight —2 density. 

The trace-free part of the Codazzi-Mainardi equation shows that the trace-free curl 
of the trace-free second fundamental form is a weight 1 conformal invariant: 

( 2 -9) Wj bc = 2V {[c n fe]a}o . 

Finally, we will need a certain novel four-manifold hypersurface invariant: 

Lemma 2.1. Let d = 3 and 

ft ft _ ft 

Then defines a symmetric, weight —1 conformal hypersurface density 

Bab € r(0 2 T*E[-l]) 


given in a choice of scale g € c by B® b 







8 


Glaros, Gover, Halbasch Sz Waldron 


Proof. This result can be established by examining the conformal transformation prop¬ 
erties of the Cotton tensor, mean curvature and the divergence of the Weyl tensor: 

C abc = n-HCatc - w abc d r d ) T , H n2 » = n-\H + r.n ), 

( 2 . 10 ) 

where Y a := For the Cotton tensor transformation we must decompose = 

TJ + h d T.h, the second term of which cancels the term in B]) fe 9 produced by the mean 
curvature. Setting d = 3 and keeping only the symmetric part of the Weyl divergence 
terms gives the final cancellation. □ 

Remark 2.2. Invariance of the density B a fe for almost Einstein structures follows from 
|Covl01 Proposition 4.3], where it was related to the Bach tensor. It also appeared in the 
construction of a tractor analog of the exterior derivative in |CLW15] , Hence we shall 
call B a fe the hypersurface Bach tensor. * 


— 72 . __— . 


= f! 


-1 


(y c W T abc + (d - 4)T c Wj ab)c + (d - 2) T c wj ab]c ) , 


2.1.1. Invariant operators. Conformally invariant operators are maps between sections 
conformally weighted tensor bundles. These play a fundamental role in conformal and 
hypersurface geometries. Quite generally we will be interested in examples of these whose 
target and domain have differing weights and possibly tensor types. There are two main 
examples that are crucial to following developments. 


Proposition 2.3. Let d > 3 and X ab be a rank 2, symmetric, weight —d = — d — 1, 
trace-free hypersurface tensor. The mapping, given in a choice of scale by 


( 2 . 11 ) 


X ab 


t \rab . 
L ab*- • — 


X7 a X b X ab + P ab X ab , d > 3, 

X a X b X ab + Pj b X ab + HII ab X ab , d = 2 , 


determines a conformally invariant operator mapping 


r(©*T£[-d]) r {®lTT[-d}). 


Proof. In dimensions d > 3, the result is standard and germane to any conformal geom¬ 
etry from the theory of invariant differential operators based on BGG sequences, see for 
example |ES97j . When d = 2 one can choose a scale, and transform the metric explicitly 
under g ab —>• Pt 2 g ab and then directly verify conformal invariance using the mean curva¬ 
ture transformation given in Equation (12.101) as well as those for the Schouten tensor and 
a weight w vector v a (which extends to higher rank tensors by linearity): 

p "? =K b -Var b + r a r b -±r 2 g ab , 

xf g v a = n w (X g b v a + (w + l)T b v a - T a v b + 5%T.v) . 


□ 


Observe that L ab is well-defined when acting on the conformal density II we record 
an identity for this required in future developments: 

Proposition 2.4. Along hypersurfaces in a four-manifold one has 


abfr2 


L ab n 


(ab)o 


= i(v a n bc )(v a ni, c ) + i(v.n a )(v.n°) + ln a6 An af) - jfl k - |n ab v c w^ bc + W abc w c 


2 fiab a 


4frabT7cwT 


rabc ti/T 
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Proof. The proof is a combination of basic tensor algebra and repeated usage of the 
trace-free Codazzi Mainardi Equation (12.9[) applied to the double divergence of 
We record the main steps below: 


v a Vn 2 afe)o 


= 2v a (n c(a v 6 n^) - 1 -Ak 


= (v.n a )(v.n a ) + 2n ab v c v a n fec + (v a n bc )(v b n ac ) + 

= ^(v a n 6c )(v a n bc ) + i(v.n a )(v.n a ) - ^n afe An ab + 
+ 2n ab v c v a n bc + n ac [v a ,v b ]nji 

= i(v a n 6c )(v a n bc ) + i(v.n a )(v.n a ) + h ab An ab + 

-ifi“ b v c wj bc + in“ c [v a ,v b ]nt 


n ac [v a ,v b ]n^- 
(v a n bc )wJ ab 

2 W Wj bc 



The third line was achieved using identity 


( 2 . 12 ) n a6 v a v.n b = | n ab An ab - |j k- 2 tr(n 2 P) - ab x c wJ ac , 

which can be proved by using Equation (12.911 to rewrite V c V a II bc = All ab — V b V.II a + 
lower order terms. In three dimensions, the Weyl tensor vanishes so the last term in 
the computation of V a V 6 II 2 afe ) 0 displayed above Equation (I2.12p can only contribute 

intrinsic Schouten and scalar curvature terms; these are precisely — P afe II 2 ; ^ 0 — |j K and 
this completes the proof. □ 


The second conformally invariant hypersurface operator needed implements the Robin 
combination of Neumann and Dirichlet boundary conditions and was shown to be con¬ 
formally invariant by Cherrier: 

Definition 2.5. Given a weight w conformal density / £ T(£M[w]), we define the 
(conformal) Robin operator 

(2.13) 6nf := V^/| s -i/;iT/| s . 

Remark 2.6. The operation (V* — yV.n)/ is the Lie derivative of the density / along 
the unit normal vector n, so by construction S^f is a weight w— 1 conformal hypersurface 
density. * 

The Robin operator can be extended to also act on conformal tensor densities, in what 
follows we will need the case of a rank two, symmetric conformal tensor density: 

Lemma 2.7. Suppose X ab £ T(© 2 T* M[w]) obeys h.X b |s = 0 then one has a conformally 
invariant operator given here in a scale g £ c 

(2.14) 5 h X ah := (V^X ab - (w - 2 )HX ab f | E . 

This defines a section of T(Q 2 T*T,[w — 1]) in the corresponding induced scale g £ c. 

Proof. This is a member of a general class of results that can be obtained from the 
Laplace-Robin operator discussed in Section 15.51 A simple proof of the statement uses 
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the conformal transformation property of the mean curvature given in Equation (12. 1011 
and that for the normal derivative of a rank two tensor density: 

vf 9 X ab = n^fVnXah + (w- 2)T ,hX ab - 2T {a n.X b) + 2n (a T.X b) ) . 

Since h.X b |s = 0, and the formula for 5^ projects onto the tangential part of XnX ab , only 
the first two terms on the right hand side above survive, and the result thus follows. □ 

3. Defining function variational calculus 

To handle hypersurface variational problems we treat hypersurfaces in terms of defining 
functions. Locally any hypersurface is the zero locus of some defining function, so no 
generality is lost studying hypersurfaces E = Z(s) for some defining function s. For 
simplicity we take M oriented with volume form cu; this and the conormal ds determine 
the orientation of E. We also shall assume, again mostly for simplicity, that E is compact. 
Hypersurface invariants are defined as follows (see [GW151 Definition 2.6]): 

Definition 3.1. A scalar Riemannian pre-invariant is a mapping 2? which assigns a 
smooth function ^(s; g ) to each Riemannian d-manifold (M, g ) and hypersurface defin¬ 
ing function s, satisfying: 

(i) ^(s;g) is natural; for any diffeomorphism (f> : M — >• M we have 

&>(<t>*s,4>*g) = t*&(s-,g). 

(ii) If s' = vs for some smooth, positive function v then 

^(•s;sOIe = ^(vs;g) | s , where E = Z(s ) . 

(iii) RZ is given by a universal polynomial expression such that, given a local coordi¬ 
nate system ( x a ) on ( M,g ), &(s,w,g) is given by a polynomial in the variables 

9abi d ai g bc i ■■■ ; d ai d a2 ■ ■ ■ d ak g bc , (detgf) )W ai ... ad , 

s, d bl s, ••• ,d bl d b2 ■ ■ ■ d bl s, \ds\~ 1 , 

for some positive integers k,l. Here 8 a means d/dx a , g ab = g(d a ,d b ), detg = 
det(# ab ) and u ai ...a d = w(<9 ai ,..., d ad ). 

A scalar Riemannian invariant of a hypersurface E is the restriction ^(E; g) := ^(s]g) |s 
of a pre-invariant &(s]g) to E. 

Remark 3.2. A hypersurface invariant & depends only on the data (M,g, E). For 
point (i) note that if E = Z(s), then 4>~ 1 ( E) is a hypersurface with defining func¬ 
tion 4>*s. In (ii), the requirement s' = vs means that s' and s are are two compatibly 
oriented defining functions with Z(s ) = Z(s' ) =: E. In (iii), since we are ultimately inter¬ 
ested in hypersurface invariants, there is no loss of generality studying defining functions 
such that \ds\~ 1 7 ^ 0 everywhere in M, since we may, if necessary replace M by a local 
neighborhood of E. 

The definition extends to tensor-valued hypersurface pre-invariants and invariants by 
considering tensor-valued functions & and requiring the coordinate components of the 
image satisfy conditions (ii), (iii), and the obvious adjustment of (i). The term hypersur¬ 
face invariant will be taken to mean either tensor or scalar valued hypersurface invariants. 
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Local Riemannian invariants, constructed in the usual way, provide trivial examples 
of hypersurface preinvariants. Examples of preinvariants that depend non-trivially on 
both s and g are given below: 


Example 3.3. Given a hypersurface defining function s, examples of preinvariants are 


&a(s\g) 


V a s 

|Vs| 


and 


^ab(s;g) = V a 



V a s V c s ^ ( V b s\ 

I Vs 11 Vs | c v |vs| y ‘ 


Upon restriction to E = 2(s), these give the unit conormal and second fundamental form 
hypersurface invariants. It is important to note that for distinguished choices of defining 
function s, these expressions simplify considerably. For example, if s obeys |Vs| = 1, 
then the second fundamental form is Jd ab = (V a Vfts)|s, but the quantity V a V;,s is, of 
course, not a preinvariant. * 


We now consider integrals dAg • over a hypersurface, where dAg is the volume 
element of the induced metric g ab . Let s be a defining function for E. Then, since (see 
Equation (12. ID ) 


choosing coordinates (x a ) 


9= [9 


ds <g> ds 
|ds| 2 


(s,y*) where g(ds,dy l )| s = 


0 , it follows that 


Thus we have the following: 


\/det i 


\j det j 

Idsl 


Proposition 3.4. Let E be a smoothly embedded, compact hypersurface in M, with defin¬ 
ing function s, and let TA be a hypersurface invariant for E with preinvariant LA(s]g). 
Then 

(3.1) [ dAgfi? = f dV g |Vs| S(s) ^(s-,g ), 

J s Jm 

where g is the pullback of the metric g to E. 

In the above S(s) is the Dirac delta distribution. The above formula is standard, see 
for example |QF03| for its use in the theory of implicit surfaces. This has the advantage 
of allowing us to perform surface integrals over hypersurface invariants & in terms of 
(often far simpler) extensions to M. 

Remark 3.5. The distributional identities xS(x) = 0 and S(f(x)) = Yl xo eZ(f) S \f{xo)\ 
for functions / with simple roots, imply |V(vs)|<5(us) = |Vs|5(s). Thus, the integrand 
\X7s\5(s)LP(s-, g) on the right hand side of (13.11) can be loosely viewed as distributional 
hypersurface preinvariant. * 

Moving now to variational considerations, let TA be a hypersurface invariant. The 
variation of the functional /(E) = f^dAg&, with respect to the embedding of the 
hypersurface E, is defined by 

= 4/ dA- gt ^ t -g) . 

0 dt y St t =o 

Here E^ denotes a smooth, one parameter family of hypersurfaces such that Eo = E 
and Ei = E outside a compactly supported subset of M. Also, gt is the pullback of the 
metric g to Ej. 
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Our strategy is, without loss of generality, to consider a family of hypersurfaces Ef = 
Z(s t ) where s t is the family of defining functions 

(3.3) s t := s + tu 


for some smooth, compactly supported Ss := u £ C°°M. Thus /(E^) = I(Z{st)). In 
general we will use S to denote the operation f(s) i—> (if/If+lhl) 

Using Equation (13.111 the variational formula (13.21) becomes 


SI = 


dt 


dV g \Vs t \d(s t )^(s t ;g) 


M 


t=0 


This reformulation leads to the following Lemma which is needed for our key variational 
Theorem 13.111 and is also a useful computational tool: 


Lemma 3.6. Let & be a hypersurface preinvariant with corresponding hypersurface in¬ 
variant PA . Consider a compactly supported variation Ss of a defining function s for a 
hypersurface E. Then 


(3.4) 


dAg PA ) = f dVg |Vs| 6 (s) SPA — — (V^ + V.n 


IM 


Vs 


where h := (Vs)/|Vs|. 


Proof. Let x be a (fixed) compactly supported cutoff function taking the value 1 on an 
open neighborhood of the support of Ss. Then 


dA- g & ) = S 


dV g x |Vs| S(s) IP 


'M 


(3.5) 


= / dv aX 


IM 


(' VfiSs)S(s ) PA + | Vs| S'(s ) Ss PA + |Vs| S(s) SPA 


= / dVgX — {V-fi)Ss 5(s) LA — 5(s) SsVnLA + \Vs\8(s) SLA 
Jm l 


Here we used that the derivative of y vanishes where Ss has support and V,^(s) = 
| Vs|<T(s). Also, here and in the above, S'(s) can be interpreted in the usual distributional 
sense. In the second line above we performed an integration by parts using that the 
integrand has compact support. The result follows using Equation (13.111 . □ 


Remark 3.7. When PA (and thus also LA) is a density of weight —d, the integral dAgLA 
is conformally invariant. In this case the operator appearing in the second term of 
Equation (13.41) is the Lie derivative of PA along h. Along E this becomes 


Vn + V.n 




where Sfi is Cherrier’s conformally invariant Robin operator; see Definition 12.51 This 
motivates the notation of the definition that follows. * 


Definition 3.8. Let E be a hypersurface and let / £ C°°M (or a density of any weight). 
Then we define 


SJ = Wnf\^-kHf, 


where / = /|s and k £ C. 
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To compute the varational gradient of dA g it remains to express the first term 
on the right hand side of Equation (13.411 in the form f M dV g |Vs|5(s) 5s PA' with some 
local formula for PA' where 


8s := 


5s 

Vs 


We have introduced the quantity 5s for homogeneity reasons: Observe that the inte¬ 
gral dAg&t being differentiated on the left hand side of (13.41) is independent of the 
defining function s t , so in particular is unchanged upon replacing s t with vs t for any 
positive, smooth function v. Thus <5( dAg&^j =: 5I(s,5s ) = 5I(vs,v5s). It is easy to 
verify the identity 


v5s 


IVMI 


5s 


|V(s)| 


so that along £, 5s enjoys the same homogeneity property as 51. To rewrite 51 with 5s 
undifferentiated, we will need the following “integration by parts” lemma: 


Lemma 3.9. Let s be a defining function and suppose C a E T(TM). Moreover let B be 
a test function (i.e., smooth with compact support). Then 


IM 


dV g \X7s\8(s)C a VjB = - 


/■ 


T /oT a 


dV g \X7s\8(s)BX7fC 


where V T = V — hV^ ; h a = (V a s)/|Vs| and C T := C — hh.C. 


Proof. First, not forgetting the dependence of V T on h, we use that B is a test function 
and integrate by parts to compute 


[ dE 9 |Vs|5(s)C a V ( [ J B = - [ dV g B 
J M JM 


Vj(|Vs| 8(s)C 1 a ) - |Vs| <5(s)C 1 b V) ( h a h b ) 


Ua 


TbwT/ 


Next, a straightforward computation shows that 


V nfl a 


vliv« 

I Vs I 


Using the distributional identity V T 5(s) = 0, the result follows. 


□ 


Remark 3.10. When C a and B , respectively, take values in sections of some tensor bundle 
over M and its dual and these are appropriately contracted to form a scalar, the natural 
extension of the above Lemma follows immediately. * 

The variation 5PA of a preinvariant is linear in the variation u = 5s and its deriva¬ 
tives up to for some finite order k. If all derivatives on the variation u are tangential 
derivatives V T n, (V T ) 2 u, ..., (V T ) fc u. we can use the above lemma to integrate these 
by parts. In fact, as we shall prove in Theorem 13.111 this is always the case. That is, 
we show there is a well-defined notion for the functional derivative of a preinvariant 
determined by 

r r A 5 c? 

(3.6) / dV g |Vs| 8(s) 5& := / dV g \ Vs| 8(s) 5s ——, 

Jm Jm os 

for compactly supported variations 5s. We now give our main variational theorem: 
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Theorem 3.11. Let £ be a hypersurface with defining function s and let PA be a hy¬ 
persurface invariant with corresponding preinvariant LA. Then for compactly supported 
variations, 

5 J dAgiiP = J dV g |Vs|<5(s) 5s - (V ft + V.h)^)) , 
where n = (Vs)/|Vs|. 

Proof. Thanks to Lemma 13.61 we need only verify the well definedness of ^ in Equa¬ 
tion (13.6p . Our first step is to construct coordinates ( s, x l ) where i = 1,..., d such that 
the vector field is orthogonal to constant s = e hypersurfaces £ £ := Z(s — e). This is 
achieved by integrating the vector fields g~ 1 (ds, •) and using the flow along their integral 
curves to push forward coordinates x % on £ = £ £=0 to £ £ . 

Since 5s has compact support, we may write the variation 5s = u € C°°M in the 
coordinates (s, x l ) as 

u(s , x l ) = u(0, x l ) + sU(s, x l ), 

for some smooth, bounded, function U. Now we decompose the one parameter family of 
defining functions in Equation (13.3p according to 

s t = s' + tu 0 , 

where the function uo is defined in these coordinates by uq = u{ 0,x*). Now note that 
1 + tU(s,x l ) is a positive function for t small, because U is bounded. (So s' = s(l + 
tU(s, x 1 )) is also a defining function.) Thus for all such t we have that /(£) := dAgLP 
satisfies 

/(£) = I(Z(s)) = I(Z(s ')). 

Now from the construction of / it follows that /(ii,it) := I(Z(s + t\sU + £ 2 ^ 0 )) is 
differentiable as a function on 1R 2 . From this and the previous observation it follows that 

dl(s + tu 0 ) 

01 =- - - 

dt t =0 

Next, from the definition of a preinvariant it follows that the variation 51 is a sum of 
terms of the form 

f dV g \Vs\d(s)C a ' - ak V ai ---V ak u, 

Jm 

for some smooth tensors C ai "' ak (s; g). However, the above coordinate argument shows 
that we may replace this by 

[ dV g |Vs|<5(s) C' ai '''“ fc V ai • • • V Ofc u 0 . 

JM 

However, since V^/ao = 0, we have Vuo = (V — hVn)uo =: V T uo- An easy induction 
shows that we can reexpress the integral as a sum of terms 

[ dV g \Vs\S(s)C a ' - a ‘Vj 1 ---Vj e u, 

Jm 

where we written u not uq because V T s = 0 and s6(s) = 0 (and C ai '" ai are some 
smooth tensors). Integrating all V T ’s by parts according to Lemma [3.91 then establishes 
Equation 13.61 □ 




Variational Calculus and the Singular Yamabe Problem 


15 


Corollary 3.12. Let £ be a hypersurface with defining function s and let 88 be a hyper¬ 
surface invariant with corresponding preinvariant 88. Then 


58 ? 

5s 


58 ? 

— 5^88 and — -(V^ + V-n)^ 2 * 


5s 

are a hypersurface invariant and its corresponding preinvariant. 

Example 3.13. Let h a = V a s/|Vs| where s is a defining function for £. Then 

8 s (s\g) = V a h a 

is a preinvariant for d times the mean curvature. Noting that 

VlJs 


5h a = 


IVsI 


it follows that 


588 = V° 


°(^p) = [A t — |V„ log |Vs|| 2 + (A t log |Vs|)] <5s, 

where Vj := V a — h a V^ and A T := 3 ab VjVj. Notice 588 does not contain normal 
derivatives acting on 5s. Thus 


5P 

— = A t log |Vs| — |VJ log |Vs 
ds 


This is not a preinvariant, but it is not difficult to verify that ^ — (V„ + ^.n)88 is and 
along £ this equals the hypersurface invariant 

— tr II 2 + d 2 H 2 + Ric(h, n) = Sc — Ric(h, n) — Sc . 

Here the last equality relied on Equation (12.3() . * 


3.1. Variational methods and identities. We now explain our strategy for computing 
variations and collect some key identities. Using Theorem 13.111 hypersurface variations 
are given in terms of preinvariants. To express preinvariants in terms of standard hyper¬ 
surface invariants, we take advantage of property (ii) of Definition 13.11 which allows us 
to choose any defining function when evaluating preinvariants along £. In a computa¬ 
tional context we are often concerned with expressions built from a sum of terms which 
separately are not preinvariants even though their sum is. Therefore in calculations it is 
important to use the same choice of defining function for every term. 

Given any hypersurface in a Riemannian manifold there is locally a defining function s 
satisfying 

|Vs| = 1. 

Defining functions with this property are dubbed unit defining functions. Since unit 
defining functions are determined by the hypersurface embedding, they provide a good 
choice for computations. In fact, there is a recursive solution to the formal problem of 
finding a smooth unit defining function si given a defining function s for a hypersurface 
such that 

I Vsi| = 1 + s e f, 

for t € Z>i arbitrarily large and / € C°°M (see |GW151 Section 2.2]). This explicitly 

relates the jets of s to hypersurface invariants. We employ the notation = for equalities 

s . 

relying on this choice while = is used when evaluating quantities along £. 

In general we abuse notation by using the same symbols for preinvariants as for their 
corresponding hypersurface invariant when the former has been explicitly chosen. Let 
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us tabulate some commonly used choices for preinvariants as well as their expressions in 
terms of a unit defining function: 



So for example, off the hypersurface the notation h a means (V a s)/| Vs| and an expression 
like X n H denotes V„Vjn a . The label T is employed to denote projection of tensors 
on M to their part perpendicular to the unit vector to, and V T := V — nV,,. 

Our first key variational identity was given in Example 13.131 but we repeat it for 
completeness: 

r-n X a 5s n a „ . i „ T . s - _ , 

(3.7) 8ri = ——j--—= Vlds = W a 5s , 

to nr 


where <5s| : = 


<Me = 


5s |. Thus 


(3.8) 


{5A ab ) T 




V> 

Ini 


T 


Vfih b 


Z V 0 V(,(5s|. 


Importantly, in the above computations, we first varied preinvariant expressions for in¬ 
variants and only thereafter specialized these to a unit defining function. 

In a functional setting, we can use Equation (13.811 to give a useful lemma for computing 
variations of functionals involving the second fundamental form: 


Lemma 3.14. Let K ab £ T(© 2 TE) and s be a unit defining function with <5s| = ds|s 
where 5s has compact support. Then 

(3.9) [ dAg K ab (dn a6 ) | s = [ dA- g 5s\V a V b K ab . 

J E Jr. 

Finally, one more technical lemma is required: 


Lemma 3.15. Let X ab be a rank two, symmetric tensor on M satisfying h a X ab = 0. 
Then 

[VfiXab - 7a&7 Cd V nXd ~ Vn(^( a fe) 0 )] T = 0. 

We leave the proof to the reader. 


3.2. Examples. As a first example, we consider the problem of extremizing the area of 
an embedded hypersurface. The relevant functional is 

L = jjAg. 

Using Theorem 13.111 we have 

51 = — f dAg <5s| (<5 - 1) = — d ■ f dAg <5s| H , 

Jt. ’ is 

which recovers the standard vanishing mean curvature H condition for minimal surfaces. 
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A second example is the variation of the Willmore energy for surfaces (in Lorentzian 
signature, this is the rigid string action of |Pol 86 | ): 

(3.10) l2 = ~JdAgK, 

where K is the rigidity density of Equation (12. 6 p . Using Theorem 13.61 we find 
sx 2 ±-±J dAg (SK) Is + \J dAg 5s\5 n K . 

Using Equation (13.91) . and observing that only the variations of Jl ab contribute when the 
preinvariant K is written in terms of H ab and 7 ab , the first of the terms above becomes 

~lfjA- g (AA)|s L- l -J^ dA - g 5s\V a V b il ab . 

The second term requires that we compute a normal derivative of K: 

5 n K - 2 HI< = V n K = 2n a6 V rt n a6 = 2n“ b [V n , V a ]n 6 
= 2R ab [ - Pa6 - - 2Era a6 ] 

= —2tr(IlP) — AHK , 

where we have used Equation (I2.13p as well as that the Weyl tensor vanishes for d < 4 
and for surfaces trll 3 = 0. Assembling the above we find 

51o = dA-g 6s\ (V a V b + Pj b + HU ab )n ab . 

Thus, using Equation (12.111) . we obtain the following: 

Proposition 3.16. With respect to hypersurface variations the gradient of the Willmore 
energy X 2 is 

--L ab R ab . 

3 

Thus the Willmore equation is L ab Jl ab = 0. Note that the Willmore invariant displayed 
in the above Proposition is manifestly conformally invariant because L ab is the invari¬ 
ant operator of Proposition 12.31 For Euclidean host spaces, using Equation (1 2. 0 p . the 
Willmore equation takes the more familiar form 

AH + 2H(H 2 — K) = 0, 

where K = ^ Sc is the Gaufi curvature. 

4. The conformal four-manifold hypersurface energy functional 

Our main variational result, as stated in Proposition 11,21 is that the obstruction den¬ 
sity ,83 for spaces agrees with the variational gradient of the rigid membrane functional 



where the membrane rigidity density L is given in Equation (12.81) . In this section we 
prove Proposition 11.21 

We begin by using Theorem 13.61 to write the variation as 

S1 3 = ^ J dA s (<5T)| S -^ j dAgds\5 h L. 








18 


Glaros, Gover, Halbasch Sz Waldron 


To compute the first of these terms, we use the Fialkow Equation (12.71) to write 

L = tr n 3 — W(h, H,n), 

so that 

^ dA- g (5L) | s = ^ dA-g [3E 2 b cffl ab - (6& ab )W ab - 21l ab W bac 5n c ] | E . 

Using Equations (I3JD and m , this becomes 

dA-gSLl ( dA- g 6s\ [3V a V 6 E 2 ab)o - W a V b W ab + 2V c (n a 6 TU^ oc )] | E 

J S 

1/ _ _ i o n — — , n n ~r n , o , _ n _i_ 

= j dA- g 5s\ [3V a V 6 nJ ab)o - V“V b LF ab + Wj ac W bac + 2E ab V c Wj ac ] | s , 


L 


where the final two terms in the last line were obtained using Equation (12.91) to rewrite 
the final term of the first line. 

To compute the second term in 51%, we employ Lemma 13.151 to first calculate the 
normal derivative of the second fundamental form: 


(VniU) T = (v n n ab )' = ([v n , v a ]n b )' = w ab - p ( ' ab)o - e^ )o 
= 21U ab - P (ab)o - 21l 2 ab)o - H\l nh , 

where in the final step we have used the Fialkow Equation (12.71) . Additionally, we must 
compute a normal derivative of the Weyl tensor: 

« i j n i j n 

V n W ab = n d n c V c W bad = ( g dc - g dc )V c W bad 
= n dc W dabc + n c V d W dabc - V J W ba c 


oj 71 _ n o n 

= n dc W dabc - 3 HW ab + n c C bca - V c lU b T ac + E c a W cb 
Orchestrating the above results we find 

1 0 0 °L ° 1 71 O 7 71 

5 n L :=V n L + 3 HL = 3E 2 b V n E a6 - (V n Jl ab )W ab - E ab V n W ab 


= 3P a 6 n? ab)o - W(n, P, n) - HW(n, E, n) - C(n, E) - V C E ab Wj ac 

- 7W(n, E 2 , n) + K 2 + R ab tl cd W dabc + 2 W ab W ab . 

Using these, the variation 61 .3 becomes 

+ 3P ab E ( ab ) 0 - V a V b W ab - W(n , P, n) 


^3 = \ [ dA s 5s\ 


3 V a v fc n 2 ab)o 


+ n ab V c Wj ac - C(n , E) - HW(n , E, n) + Wl ac W 


n n 


bac 


- 7W (n, E 2 , n) + K 2 + E ab E cd W dabc + 2 W ab W ab . 

This result can be written in a manifestly conformally invariant way using the opera¬ 
tor L ab of Proposition 12.111 and the hypersurface Bach tensor of Lemma 12.11 

1 P ^ / 0 ™ o 7 _ 71 -1- 71 

5 X 3 = - ^ dA-g 5s I [L ab [3E 2 a6)o - W ab ] - E afe B ab + K 2 + W T bac W bac 

- 7W (h, E 2 , h) + II“ b E cd TU dabc + 21U ab lU afe ) . 
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This shows that the tensor 

B 3 = - . 

■’ 6 L 


L ab (ml h) -w ab ) - n ah B ah + K 2 - 7 w ab ul h + 2 w ah w ab + ir h n cd w cabd + w^w 

is the variational gradient and completes the proof of Proposition 11.21 


abc 


5. Conformal hypersurfaces and the singular Yamabe problem 


We define here a singular Yamabe problem as the boundary version of the classical 
Yamabe problem of finding a conformally related metric with constant scalar curva¬ 


ture [AM881 IMaz91l IACF92] : 


Problem 5.1. Given a d > 3-dimensional Riemannian manifold ( M,g ) with bound¬ 
ary £ := dM, find a smooth real-valued function u on M satisfying the conditions: 

(1) u is a defining function for £ (i.e., £ is the zero set of u, and du x 0 Vx € £); 

(2) g° := u~ 2 g has scalar curvature Sc 5 = —d(d — 1 ). 


Where 0 < u := p 2 ^ d 2 ), part (2) of this problem is governed by the Yamabe equation 


_ 4 ^1 A » + Sc» 


d +2 


p + d (d — 1) p d - 2 = 0 . 


But since u is a defining function, to deal with boundary aspects the equation can be 
usefully recast: Since each g E c is in 1:1 correspondence with a true scale r £ T(£Af[l]), 
setting the smooth defining function u = a /r, the Yamabe equation becomes 

(5.i) + 

In the above A := g ab X ’ a V& and all index contractions are performed with the conformal 
metric g and its inverse. It follows easily that the quantity S(a) on the left hand side of 

Equation (15.111 is conformally invariant. Note that S(a) = — Also, because u is 

a defining function, a is a defining density : this is a section of T{£M[ 1]) with zero locus 
2(a) = £ and such that for any Levi-Civita connection in the conformal class Vcr 0 
along £; see Section T5.3I below. 

Development of the asympototic analysis of the singular Yamabe problem as well as 
uncovering its fundamental role in the theory of conformal hypersurface invariants uses 
the tractor calculus treatment of conformal geometries. We next review key aspects of 
that theory. 


5.1. Conformal geometry and tractor calculus. Although there is no distinguished 
connection on the tangent bundle TM for a d-dimensional conformal manifold (M, c), 
there is a canonical metric h and linear connection V^" (preserving h) on a related higher 
rank vector bundle known as the tractor bundle TM. This is a rank d + 2 vector bundle 
equipped with a canonical tractor connection V 7 " |Tho26l [BEG94 ], The bundle TM is 
not irreducible but has a composition series given by the semi-direct sum 

T a M = £M[l]T£ a M[l}T£M[-l ]. 

There is a canonical bundle inclusion £M[— 1] —>• T A M , with T*M[1] a subbundle of the 
quotient by this, as well as a surjective bundle map TM —» £M[ 1]. We denote by X A 
the canonical section of T A M[ 1] := T A M <S> £M[ 1] giving the first of these: 

(5.2) W 4 : SM[- 1] -> T A M , 

and term X the canonical tractor. 
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A choice of metric g € c determines an isomorphism (or splitting ) 


(5.3) 


TM k £M[ 1 ] © T*M[ 1 ] © £M[- 1 ] . 


We may write, for example, U = (a, p a , p), or alternatively [ U A ] g = (a, p a , p ), to 
mean that U is an invariant section of TM and (a, p a , p) is its image under the 
isomorphism (15.311 : often the dependence on g € c is supressed when this is clear by 
context. Changing to a conformally related metric g = Q?g gives a different isomorphism, 
which is related to the previous one by the transformation formula 

(<r, Pa,pf 9 = {Va, £l(p a + <rT a ), n~\p - g bc T b p c - \ag bc T b T c ))\ 


where T 5 is the one-form 1 dfl. 

In terms of the above splitting, the tractor connection is given by 




( a \ 

( V a (7 - Pa 

\ 

(5.4) 

vr 

Pb ■■= 

Tapb T g a bP T Pab& 




{p) 

\ V a p - P ac p c 

/ 


We will often recycle the Levi-Civita notation to also denote the tractor connection; it 
is the only connection we shall use on TM , its dual, and tensor powers. Its curvature is 
the tractor-endomorphism valued two form 77^; in the above splitting this acts as 


(5.5) 


/ a N 


( (T ^ 


0 ) 

Pc 

V p ) 

= [v a ,v b ] 

Pc 

V p ) 

— 

fT abc pd + C abc <7 
\ ~CabcP C ) 


For [ U A ] = (cr,p a ,p) and [ V A ] = (r, v a ,K,), the conformally invariant, signature (d + 
1,1), tractor metric h on TM given by 


( 5 . 6 ) h(U, V ) = h AB U A V B = ctk + g ab p a v b + pr =: U-V , 

is preserved by the tractor connection, he., V^/i = 0. It follows from this formula 
that X A = 1iabX b provides the surjection by contraction i(X) : TM —>• £[!}■ The 
tractor metric h A B an< 3 its inverse h AB are used to identify TM with its dual and to 
raise and lower tractor indices. The notation V 2 is shorthand for V A V A = h(V,V). 

Tensor powers of the standard tractor bundle TM , and tensor products thereof, are 
vector bundles that are also termed tractor bundles. We shall denote a tractor bundle 
of arbitrary tensor type by T'^M and write T^M[w] to mean T®M ®£M[w]; w is then 
said to be the weight of T®M[w\. 

Closely linked to V 7 " is an important second order conformally invariant differential 
operator 

D A : T(T®M[w]) -> T(T^M[w - 1 ]), 

known as the Thomas (or tractor) D-operator. Here T®'M[w — 1] := T A © T®M[w — 1]. 
In a scale g, 


/ (d + 2 w — 2 ) w ^ 

(5.7) [D A ] g =\ {d + 2w — 2)V a , 

V -(A9 + J w) 


where here A = g ab X a X b , and V is the coupled Levi-Civita-tractor connection |BEG94l 
lTho26j . The following variant of the Thomas D-operator is also useful. 













Variational Calculus and the Singular Yamabe Problem 


21 


Definition 5.2. Suppose that w / 1 — The operator 

D A : r(T $ M[w]) —> r(T # 'M[w; - 1]) 

is defined by 

D A T := ---- D A T . 

d + 2w — 2 

Remark 5.3. We term the critical value w = 1 — ^ the Yamabe weight, since one then has 

D A = -X A Uy , 

where the operator — Dy : T(T $ M[ 1 — |]) —> T(T®M[— 1 — |]) is the conformally 
invariant, tractor-coupled, Yamabe operator given by A — J in a choice of scale. * 
The Thomas D-operator is null, in the sense that 

D a o D a = 0; 


this operator identity can easily be derived from Equation (15. 7 j) . The Thomas D-operator 
is not, however, a derivation. Its failure to obey a Leibniz rule is captured by the following 
Proposition, a result which was first observed in [ JTW13] and proved in [GW15] : 

Proposition 5.4. Let T t £ T(T®M[w{]) for i = 1,2, and hi := d + 2 Wi, h \2 := d + 
2w\ + 2w2 — 2 with hi 0 ^ h\ 2 . Then 

(5.8) D a (T ] T 2 ) - (D a T 1 )T 2 - T!(D a T 2 ) = -- Y a (. D b T 1 ){D b T 2 ). 

d + 2wi + 2 w 2 — 2 

We shall also need the following Lemma which is easily verified by direct application 
of Equation (15.71) or the modified Leibniz rule (I5.8|> in tandem with the identity 

(5.9) D A X b = h AB . 


Lemma 5.5. Let T € T(T $ M[w\). Then 


(5.10) 


D a {X a T) = {d + w)(d + 2w + 2)T . 


5.2. Conformal hypersurfaces and tractors. The basis for a conformal hypersurface 
calculus [GovlOl IGra031 ISta051 |Vyal3[ IGW15| is a unit tractor object N A £ T(TM)|s 
(from jBEG94j l corresponding to the conformal unit conormal h, termed the normal 
tractor. This is defined along £, in a choice of scale, by 


(5.11) 


[N a ] 


( 0 \ 

\ h a N a N a = 1. 

\- H J 


The subbundle TM t orthogonal to the normal tractor (with respect to the tractor met¬ 
ric h) along E is canonically isomorphic to the intrinsic hypersurface tractor bundle TY, 
see |BG01| and |Govl0l Section 4.1]. This is the conformal tractor analog of the Rie- 
mannian isomorphism between the intrinsic tangent bundle TE and the subbundle TM t 
of TM |y orthogonal to h a along E. This isomorphism identifies these bundles and thus 
we use the same abstract index for TM and T E. In explicit computations in a given 
choice of scale, we will need to relate sections given in corresponding splittings of the 
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ambient and hypersurface tractor bundles. In terms of sections expressed in a scale g E c 
(which determines g E c), the isomorphism is given by 


(5.12) 

M, : = 

/v+\ 

V a 

* 

( v+ \ 
v a - n a Hv + 



\V~J 


\y~ + \H 2 v + j 





where E" 4 E T(TM t ), E" 4 E r(T£) and the SO(d + 1, l)-valued matrix 



( 

1 

0 

°\ 


h a H 


0 

V- 

-W 2 

h b H 

1/ 


In the above the canonical isomorphism between T*M t and T*£ defined by the unit 
normal vector h is used to identify sections of these bundles. 

The ambient tractor connection V obeys an analog of the classical Gaufi formula 
relating it to the intrinsic tractor connection V of (£,c) which is known as the Fialkow- 
Gaufi formula ( |Gra03l !Sta05l Vyal3| or in our current notation |GW15l Proposition 
7.14]) because the Fialkow tensor of Equation (12.711 encodes the difference between these 
two connections along £. As in the Gaufi case, the tractor-coupled operator V T : = 
V — nVft is well-defined along £ irrespective of how tensors are extended to the ambient. 
For future use, we record some identities for this operator acting on the normal tractor: 


Lemma 5.6. Along £, tangential, tractor-coupled gradients of the normal tractor are 
given in a choice of scale by 




( 

* \ 


( 

-n ab 

\ 

(5.13) 

V T a N c = 


n ac 

, vTv 4 'iv c = 


Ilf 1C — j_^5acV.Ilf, 




\~ 

v.n„ 



1 V7 TT, TT c P t 

/ 



d-1 / 




while the tangential, tractor-coupled operator A T := < 7 a6 VjVj~ gives 


(5.14) 


A t N c 


( 0 \ 

fffV.H C - n c K 

v.v.n+(d-i)p at ,n a6 
\ d-1 / 


Proof. Since V T is tangential, the first identity follows directly from |GW151 Corollary 
6.7]. The third identity follows directly from the second, both of which only require a 
direct application of Equation (15.4p . □ 

Remark 5.7. Note that expressions such as Vjllbc in the above Lemma are well-defined 
along £ since V T is tangential. It is obviously possible to further develop the second 
and third identities using the Gauh formula (12.21) and (for d > 3) the Fialkow-Gau£ 
Equation (12.71) . * 


5.3. Singular Yamabe problem asymptotics and defining densities. Given a hy¬ 
persurface £, a section a E r(£[l]) is said to be a defining density for £ if £ = Z(a) 
and Vo - is nowhere vanishing along £, where V is the Levi-Civita connection for some 
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g E e. The conformal analog of the exterior derivative of a defining function is then 
provided by the tractor field 


(5.15) 


la '■= D A a = (a, V a cr, -^(A + J)cr) =: (a, n a , p). 


In Riemannian signature for any defining density a, we have that 

Il> o 

holds in a neighbourhood of £. More generally, any section a € T(fM[l]) such that I A 
is nowhere vanishing is called a scale, I A is its scale tractor and the data (M, c, £) 
is an almost Riemannian structure (or geometry). The following proposition, whose 
proof follows directly from Equations (15.151) and (15.61) . shows that almost Riemannian 
geometries are a natural setting for the singular Yamabe problem: 


Proposition 5.8 (see [Gov TO]). For a € T(£[l]) the quantity S(a) of Equation \5.1\) is 
the squared length of the corresponding scale tractor: 

(5.16) S(a) = II := h AB I A I 
Indeed, the Yamabe equation (15.11) now reads simply 

Il = I- 

Thus, the asymptotic version of the singular Yamabe Problem 15.11 is stated as follows: 
Problem . 5.9. Find a smooth defining density a such that 

(5.17) ll = \ + a t A t , 

for some smooth Ag € T (£[—£]), where £ € IN U oo is as high as possible. 


The above problem was solved in [GW131 IGW15| . that result is captured by the 
following Theorem whose proof may be found in [GW15] : 

Theorem 5.10. Let a be a defining density for a hypersurface £ and a q := a/^/7%, then 
there exist smooth densities Ak € T(££[— k]) such that 

(5. 18 ) =■- 1 1 1 71 " 1 71 1 1 71 

solves 


O — (Jq ^1 + AiCq + A 2 CTQ + • • • + 


= 1 + ctq B a , 

for some smooth B a € T{£M [—d]). For 1 < k < d the densities are determined by 
the recursion 

d H . ~ 1 


ffc — Gk -1 


1 - 


CT fc- 1 


2 (d-jfe)(lfe + l) 


where (Jk-l '■= ^(l + Y^i a ^j) solves 7| fc _ 1 = l + <r fc 1 C for smooth C € T(£M[1 — k]). 
Moreover, for any true scale r E T{£M\ 1]), 


a = a 


1 + 


solves 


la' = 1 + ° d+l (C r + D log(cr/r) + D 1 log 2 (ir/r)) 


for smooth C',D,D’ € T{£M\—d — 1]). Finally, given two defining functions a and a 
for £, then 


= Ba\ y =: Bd . 
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The last statement of this theorem implies that the density is uniquely determined 
by (M, c, X). Since it obstructs a smooth all orders solution to the singular Yamabe 
problem we term it the obstruction density. The solution cr(er) is determined uniquely 
by the recursion up to the order of the obstruction. Thus, up to the given order, a is a 
canonical scale defining the hypersurface X and whose scale tractor has unit length to 
accuracy a d , so is termed a conformal unit defining density. 


5.4. Conformal hypersurface invariants and canonical extensions. A Rieman- 
nian hypersurface invariant <^(X;g) := &(s;g)\z with the property <5P(f ls,fl 2 g) = fl w 
^(s; g) is termed a conformal hypersurface covariant. This determines an invariant sec¬ 
tion a,g ) G T(£M[w}) where a is a defining density. Thus ^(a, g)|s called a 
conform.al hypersurface invariant. As for the Riemannian case, this definition extends to 
tensor and therefore also tractor-valued invariants. 

Since the conformal unit scale a is a unique defining density, up to the addition of 
terms a d+1 S for S G T(£M[—d]), it follows that any conformal density g), 

where d'(cr) is determined in terms a given defining density a via Equation (15.1811 of The¬ 
orem [5T0l gives a conformal hypersurface invariant gf)|s so long as the formula 

for & involves jets of a of order < d. (In fact, is a density-valued prein¬ 

variant.) This construction provides holographic formulae for conformal hypersurface 
invariants, see |GW151 Section 6.2]. 


Example 5.11. The scale tractor for the conformal unit defining density gives a holo¬ 
graphic formula for the normal tractor 

l£\v = D A a\z = N A . 

This result was first Droved in (Gov IS] * 

Moreover, in the above example we may view the conformal unit defining density as 
providing a canonical extension I A of the normal tractor N A . A short computation (in 
a choice of scale) shows that 


Y a I (J — j V ajlh T P ab& T < tabP 


where p := — ^ (Act + Jcr) as in Equation (15.151) . Comparing this with Equation (15.1311 
gives a holographic formula for the trace-free second fundamental form 

(V a ?r fe + P ab cr + g ab p)\z, ■ 


Hence we have the following canonical extensions of the trace-free second fundamental 
form and rigidity density 

(5.19) fl afe := V a n b + Pabd + g a bP and K := (X7 a n b + P ab a +g ab p)(X7 a n b + P ab d + g ab p). 


The above formulae are now well-defined in M rather than only along X. These extensions 
are canonical up to a finite order determined by Theorem 15.101 


5.5. The Laplace—Robin operator and tangential operators. By contracting the 
scale tractor and Thomas D-operator, we can build a new operator that plays two roles: 
that of the ambient Laplace operator and of a conformally invariant boundary Robin-type 
Dirichlet-plus-Neumann operator ]Gov07| . Thus we call the operator 

I D : T(X*M[w\) —> T (T*M[w - 1]), 
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the Laplace-Robin operator. Calculated in some scale g£c 

in Qjiii 

(5.20) I-D = —a/S. + (d + 2w — 2 ) [V n - -(Act)] -— (d + w — l)crj , 

where n = Vex. Hence this is a Laplace-type operator which is degenerate along £ = 
Z(a). Indeed, along £ the Laplace-Robin operator becomes first order. In particular, 
for a conformal unit defining scale, along £ and in a choice of scale we have 


Ia-D = (d + 2w — 2) (V n — wH^j . 

Hence, for w ^ 1 — if and in the spirit of the Section 15.41 along £ the Laplace-Robin 
operator is a operator-valued holographic formula for the Robin operator 5n as given 
in Equation (|2.13p (recall that for a conformal unit defining density, along £ one has 
Vcr = h). Given a natural density, tensor or tractor field / along £, it is often the case 
that we can fix an extension / off £, that is canonical to some order using the conformal 
unit defining density, such as those given in Equation (15.191) 1. In that case we will denote 
the Robin operator by 


Sr I = 


IaDf | s , w + 1-f, 

Vn/| s + ^/, W = 1 -f. 


In light of this, we will often use the same notation for / and /. 

Combining the conformal unit defining density and the Laplace-Robin operator allows 
successive normal derivatives of canonically extended conformal hypersurface invariants 
to be defined up to the order of the obstruction. The canonical extensions of the normal 
tractor and intrinsic canonical tractor are the scale tractor I A and the ambient canonical 
tractor X A , respectively. The following Lemma gives one normal derivative acting on 
these: 


Lemma 5.12. The Robin operator acting on the canonically extended canonical and 
normal tractors gives: 


(5.21) 


d~R,X A = N a , 


SrN a 


kx a 

d- 1 


Proof. The first identity is a direct consequence of Equation (15.91) . For the second we 
compute in a scale 

( 0 \ 

6 r N a = V n I A | s = V n n a - Hn a . 


\X n p- P(n,n)y 

To complete the proof we need the normal derivatives of the canonical extensions of n a 
and p determined by I A . These were computed in | GW15 . Lemmas 6.1 and 6 . 6 ] and, 
along £ are: 


(5.22) 


^ n^a — Hn a , 


S7 nP 


P(n, n) + 


K 

d- 1 ' 


□ 


The Laplace-Robin operator has one further crucial role to play. In [GW14| . it was 
shown that by viewing the scale a as an operator mapping F(T®M[iv]) —» T(T^M[w + 
1]) and acting by multiplication, then the commutator of this with the Laplace-Robin 
operator acting on weight w tractors obeys 

[I c -D,a] = -ll(d + 2w). 
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This identity underlies a solution-generating si(2) algebra. A particular consequence of 
this is that the operator 

/I \ k 

n ■■= {-pjvD) 

is tangential when acting on weight w = tractors. Specializing to the conformal 

unit defining density gives a differential operator 


P/. : r(r*M[h-l±i]) 




which is completely determined by the data (M, c, E) |GW151 Section 8.1]. We call this 

__ k 

the extrinsic conformal Laplace operator. For k even the operator P& is (A ) 2 + L.O.T., 
up to multiplication by a non-vanishing constant, where “L.O.T.” stands for terms of 
lower derivative order involving both intrinsic and extrinsic curvature quantities. 


6. Singular Yamabe obstruction densities 

Given a unit conformal defining density, the obstruction density B f j has a simple holo¬ 
graphic formula whose leading structure is given in terms of the extrinsic conformal 
Laplace operator |GW15l Theorem 8.11]: 

Theorem 6.1. The obstruction density is given by 
2 


( 6 . 1 ) 


B, = 


dl(d + 1)! 


Da 


J B 


P 3 N B + (-1 ) d - 1 [LD d ~ 1 (X B K)] |, 


where the rigidity density K is canonically extended off E by the formula K = PabP AB 
with P AB := D A I B . Also, the projector E^ := 5^ — N a Nb- 

The main aim of this section is to compute £>3 using the above holographic formula. 
The canonical extension of the rigidity density K stated in the theorem is the same as 
that given in Equation (I5.19() . Also, in hypersurface dimensions d = 2,3, the extrinsic 
conformal Laplacians acting on weight zero tractors are given by 

P 2 = A T , 

( 6 . 2 ) 

3 0 ■ LL v a v b 


P 3 = -8 n a6 V J- 8 V.n fc V,T + 4n n (Tzf ib o V 6 + V b O n 


Proposition 6.2. If d = 2, then 

P 2 N c = -X u L ab H a0 - N U K , d = 2. 
Proof. Using Equations (16.2[f and (15.14ft . we have for d = 2, 

( 0 

p 2 N c = A t N c = —n c K 

V-V.V.H - P ab W lb 

Using Equation (12.111) the quoted result follows. 


□ 


Remark 6.3. Our computation of £>3 relies on a d = 3 version of the above Proposition. 


The second term in Equation (16.11) involves normal derivatives of the canonically ex¬ 
tended rigidity density K. The following proposition provides the first of these. 
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Proposition 6.4. The canonically extended trace-free second fundamental form.h, rigid¬ 
ity density K and membrane rigidity density L obey 

1 

T 


(6.3) 


— — n ab + Wab + — Q a bK , 

S r K = —2(J — 2)L . 

Proof. We begin by computing a normal derivative of the canonical extension of IR;,: 

^ n^-ab —V a^b T gab' V nP T nPab 4“ Pab^ 

n 

—4?ab a^c)V IT'b V fl Vb(P®0 + g a b V+ Pab T O’ (V^ 2p) P a fo 

71 o 2 

— P‘uh IRb T pHab 27i^ a Vb)P d - gab^nP T Pab 


(6.4) 


d- cr 


(v n - 3p)P ab + 2n c(a Pg } - V a V b p + 0(a 2 


Projecting onto the tangential piece of this quantity along £, and using Equation (|5.22l) 
we find 


1 


(v n n a6 ) = -n 2 ab + w ab + —g ab K - Hn ab . 


Using Equation (12.141) . we have 

(<5 R n afe ) T = [(V n + H) n 


ab | — — Ilab + W a b + 7;9ab& ■ 


The second equation follows by noting that 5rK = 2II ab <5 R n a {, and using the above 
result. □ 


Remark 6.5. Recalling that the Willmore energy functional (13.101) for embedded surfaces 
is the integral of the rigidity density K and its analog (11.ID for embedded spaces is the 
integral of the membrane rigidity density L, we see from Equation 16.31 that the Robin 
operator relates the two integrands for these functionals. It would be interesting to 
investigate whether this phenomenon holds in higher dimensions d = 2k to d = 2k + 1 

(k € Z> 2 ). * 

Using Propositions 16.31 and 16.21 as well as Equations (15.2111 and (15.10D . we arrive at an 
explicit formula for the obstruction density for surfaces. 

Proposition 6.6. The obstruction density for d = 2 is given by 

B 2 = - l -L ab il ab . 

Comparing the above with Proposition 13.161 verifies again that the obstruction density 
and the gradient of the Willmore energy X 2 agree. 


6.1. Four-manifold obstruction density. Our goal is now to find hypersurface for¬ 
mula for the d = 3 obstruction density and thus prove Proposition 11.11 We begin by 
calculating P^N B . 

Proposition 6.7. When d = 3, 

P 3 N b = AB'X b + - [U-^c^K - A(S r K)N b , 

3 
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where 
(6.5) 

B' =L ab {2t 2 (ab)o - W ab ) - ll ab B ab + -K 2 - 3W (n, II 2 , h) + W ab W ab + \w T abc W 


abc 


Proof. Using Equation (15.131) . the first two terms in Ps-ZV^' appearing in Equation (16.21) 
can be written explicitly as 


/ 


8K 


\ 


-8n a6 V a II 6c - 4H(!V.n fe + 8 n c (tr n 3 + HK) 
4H" 6 v a v.ii fe + 4v.n a v.n a + 8 tr(n 2 p T ) / 


To compute the remaining two terms of PsIV^, we make use of the tractor curvature and 
connection formulas in Equations (15.51) and (15.41) , respectively, as well as Equation (I5.13j) , 
and find 


n a V b onl b N c = 


n a ni o V bT N c = 


0 

-w abc n ab 

_ n 

,-\7 b C b -W(n,P T ,n), 

- 0 ., 

~W abc E ab 


—C(n, n) 

In the first equation we have made use of the four-manifold identity ^J b W baa i = C c d a - 
Orchestrating, we have 


p.uv G = 


( 8 K 

- 8 II ab V a f[ bc - 4n|!V.n 6 - 8 W(n, D, c) T + 8 n c ( tr II 3 - W(n, H, n) + HK) 

o _ _ o — O _ O O O _ n 

V4n a 6 V a V.n 6 +4V.n a V.n a +8tr(n 2 P T )-4C(n,n)-4V 6 C 6 -4VE(n,P T ,n) 


This expression can be simplified by comparing it with the boundary Thomas D-operator 
acting on the rigidity density, 

( 6 . 6 ) 


■P 


- 11 B 
C 


D C K = 


/ 1 

0 

°\ 


/ 8 K \ 


( 8K \ 

n b H 

\-h H 2 

-n c H 

0 

V 


-4 V C K 

{-§{&-21)K) 

H 

—4S7 b K + 8 n b HK 

|(A— 2~])K—AH 2 K J 


Here we have used the isomorphism (15.121) between sections of TS and TM t . To match 
this expression with terms appearing in P 3 N c as displayed above, we make use of the 
trace-free Mainardi Equation (12.91) to obtain the identity 

\vcK = n ab v c n afe = n ab v a n bc + ^ff^v.n, + w T abc k ab . 

Using this identity, Equation (|2.12D . the Fialkow Equation (12.71) as well as Proposition [673] 
we have 

P 3 IV- b = - [U- l } b c D c K - 4 N b SrK + 4 B'X b , 

3 
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with B' given by 

2 - - - - 4 - 4- 8 _ n 

b' = ^v a n fec v a n fec + (v.n a )(v.n a ) + -n flfc An afe - -m - -i i ab v c wj bc - 3 W(n,n 2 , n) 

+ \k 2 + W ab W ab + 2i ab \7 c wj ac - V b C b - W(n, P, n) - C(n, ft) - HW(n,TL,n ). 

It remains only to show that this last expression reduces to that stated in the proposition. 
This follows directly from Proposition 12.41 and the identities 

c b = v c it c6 + wj cb n ac => w b c b = v a v b w ab + n ab v c wj bc + -wl bc w abc ■ 

these are direct consequences of Equation m and m • Noting that n = h along S 
completes the proof. □ 

We now turn to the second term in the holographic formula (16.11) . We first need the 
following technical result: 

Lemma 6.8. Let T € 2]) and d = 3. Then 

I-D 2 (X c T) = -2 D C T. 

Proof. The key to the proof is to note that for d = 4, the tractor X (: T has ambient 
Yamabe weight w = — 1, so that 

I-D{X b T) = -aU Y {X B T ). 

Acting with the remaining I D = —25r (at this weight) and calculating explicitly using 
the tractor connection (15.4p gives the identity stated. □ 


Our result for two canonical normal derivatives of the rigidity density multiplied by 
the canonical tractor is given by the following: 


Proposition 6.9. Let d = 3. Then, 

s 


(6.7) 

where 

( 6 . 8 ) 


I-D 2 (X c I<) = AX b B" - - [U~ 1 ] B cD c K + m B 5 R K , 

3 


B" = L ab il 2 ab)o + h< 2 - 4W(h, n 2 , n) + W ab W ab + J4 ab t cd W cabd + \w T abc W abc . 


Proof. Using Lemma [6.81 we must compare 


-8K \ 


-2D C K = 




4 V C K 
\2(A - 24)K) 


with Equation (16.61) . This yields Equation (16.71) with 

(6.9) B" = ]x7 2 n K + 2 H5 R K + ^ AK - \h 2 K - i] K - ^K 2 - P(n, n)K. 


It remains to show that this expression reduces to Equation (16. 8p . For that, we need to 
compute two normal derivatives of the canonically extended rigidity density K. Return¬ 
ing to Equation (16.41) . we see that the following normal derivatives of curvature terms 
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appear in i\7 2 iP: 

n 

V n R a b + 2V n P ab 

= {V n n d )R dabe n e + n c n d V c (R da b e n e ) + 2V„P ab 
= HR ab + (, g cd - g cd )V c (Rdaben e ) + 2V n P ab 

71 7 -T ™ 7 

= HR ab + V ( Relabel ' C ) + i? b ' a + 2V n P ab 


— HRab + (V dnC )Rdabc ~ n C (V b .R c ; ac <1 + V cRdc^b) + Vj -R b d a + 2V n P, 


Tn d 


ab 


= HR ab + (t dc + Hg dc )R dabc + n c V b R ac - g ab V n J + VjR b d a . 

In the fourth line above we used the Bianchi identity for the Riemann tensor. Upon 
contraction with II ab the above yields 

tl ab {y n R ab + 2V n P ab ) = W (n, II 2 , n) - HW(n, H, n) - 3H tr(IIP) + tl ab Jl cd W cabd 

- tr(E 2 P) + fl ab V a p{ + KP(n,n) - t ab V d Wj ad . 
Combining the above expression and Equation (16.411 . we have 

n ab v 2 n ab I n ab (v n R ab + 2V n p ab ) - 2n a6 n cb v n n^ + kv hP 

+ ptl ab V n E ab - 3ptr(IlP T ) + 2tr(n 2 P T ) - n ab v a v bP 


( 6 . 10 ) 


= ^n afe An afe - \u< + 2 ETtrn 3 + \h 2 K + --K 2 - 2IU(?r,n 2 ,n) 

o o Z 4 

/I 

- 2HW(n, n, n) + fl ab JL cd W cabd + KP(n, n ) - -ir b V c IU 6 T ac . 


Here we have used Equations (15.2211 . (16.411 and (12.711 as well as the 3x3 matrix iden¬ 
tity trll 4 = t \K 2 . In addition, to handle the term tl ab 'V a V b p, the following consequence 
of the Mainardi Equation (12.41) is needed: 

n ab v a v b H = hi ab AE ab - n ab v a p b T - ^ m - tr(n 2 p) - ln ah v c wl lc . 

The computation of V 2 AT also requires the cross term (V n il ab )(V n II ab ). Again using 
Equation (16.411 we have 


o £ " „ - o o o 2 1 

T H^ab -^a6 "P ~^9abR ■ 

Contracting this identity on itself we find 

O O , 71 71 O O 0 1—0—0, 

(V n n ab )(V n n ab ) = W ab W ab - 2HW(n,P,n) - 2W(n, II 2 , n) + -(V.n b )(V.n b ) 

+ H 2 I< + 2 H tr(n 3 ) + i K 2 . 

Orchestrating this and Equation (16.101) . we have 

\v 2 n K = i(V.n o )(V.n a ) + ^II ab An ab + ^ K 2 - 41U(n,n 2 ,n) + w ab w ab + \h 2 k 

- AHW(n, ft, n) + 47f trH 3 - hK + n“ b n cd IU caM + P (n,n)K - Ul ab V c Wj ac . 

o o 
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Substituting this result into Equation (16.91) and simplifying gives 

1 _ 07_0 1—0 —O 9 O ,_0 2— 4 O 1 — n -|- 

b" = -v a n fec v a n fec + ~(v.n a )(v.n a ) + -n ab An ab - -m - -n ab v c wJ ac 


- 4 W{n, II 2 , n) + X -K 2 + W ab W ab + t ab U cd W cabd . 


Using Proposition 12.41 gives the result (16.81) . which completes the proof. 


□ 


We now complete the proof of Proposition 11.11 giving the obstruction density for four- 
manifolds: 


Proof of Provosition \l.l[ Using Equation (16.11) and Propositions 16. 7l and 16.91 the obstruc¬ 
tion density is given by 

Bz = 4d a I±X a {B' + B")}. 

Substituting the formulas (16.5p and (16.81) for B' and B" into the above and using 
Lemma 15.101 gives the result. □ 

6.2. Examples and umbilicity. Our results can be easily applied to explicit metrics 
and hypersurfaces. As we shall explain below, this gives both an independent check of 
our computations and also generates interesting examples. 

For surfaces, total umbilicity suffices to ensure vanishing of the obstruction density 
(this follows by inspection of the formula for £>2 given in Proposition 16.6p . The following 
proposition shows this is not true in dimension d = 3. 

Proposition 6.10. Total umbilicity (i.e., vanishing of II everywhere) is not a sufficient 
condition for the obstruction density to vanish. 

Proof. See the explicit counterexample on the third line of the table below. □ 

Remark 6.11. There is no in principle difficulty constructing examples of umbilic hy¬ 
persurfaces with non-vanishing obstruction density in dimensions d > 3. Also, albeit in 
Lorentzian signature, the Schwarzschild example below provides another d = 3 example 
of this. % 

In the following table we give some sample four manifold metrics, hypersurfaces and 
their obstruction density £> 3 . Additionally, d denotes the conformal unit density in the 
scale determined by the quoted metric, H the mean curvature, II the trace-free second 
fundamental form and L the membrane rigidity density. Non-vanishing expressions that 
are too long to be displayed are denoted by a *. These results were generated with the 
aid of a computer software package (PMOlj . 


ds 2 

E 

a 

H 

n 

£3 

L 

dx^ + (l-\-atx)(dy‘2+dz‘^-\-dw^) 

a;=0 

x + fx 2 -^-x 3 + ^-x 4 

a 

2 

0 

0 

0 

d x ^ + ( 1+ a x ) d y ^ ( 1 — q: x ) d z ^+ dw ^ 

x=0 

x-^-x 3 

0 

(,dy^ —dz^ 1 ) 

a 4 

U 

0 

dx“^ -\-(l-\-z)dy^-\-(l-\-w)dz^ -\-(l-\-y)dw“^ 

x=0 

~k 

0 

0 

~k 

0 

- ( l-2M) dt 2+ dP M +r 2 dn 2 

1 r 

r=R+2M 

~k 

Eq. (Ki.IIll 

Eq. CH 

Eq. (Ki.ldl 

Eq. (l().12t 


The first line of the table gives a metric and an umbilic but non-minimal hypersurface for 
which the obstruction density vanishes. In the second line, the hypersurface is minimal 
but non-umbilic. The third line is particularly interesting, the hypersurface is both 
umbilic and minimal, yet the obstruction density does not vanish. Finally in the fourth 
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line, we use the fact that all local formulas presented in this paper can be extended to 
Lorentzian structures (provided the conormal is nowhere null), and study a spherical 
shell of radius R in the Schwarzschild metric. The conormal dr is spacelike for R > 0, 
timelike for R < 0 and lightlike at the horizon R = 0, so that strictly the formulas 
presented below apply only to the spacelike case. In that case the second fundamental 
form is given by 


(6.11) H 


1 M + 2R 
3 v 'R(R + 2M) 3 ’ 




2 Rdt 2 




so that the hypersurface £ is neither minimal nor umbilic, save for the distinguished 
radius R = M where umbilicity holds. At that value the membrane rigidity density also 
vanishes; at general R it is given by 


( 6 . 12 ) 


2(Af - R){M 2 + 6MR + R 2 ) 

R^R{R + 2 M) 9 


The obstruction density is 


(6.13) 


Bs 


2 M 4 + 4 M 3 R + 30 M 2 R 2 - 8 MR 3 - R 4 
27R 2 {R + 2Mf 


and has one simple real root in the region R > —2 M at R/M ~ 2.9. 

Finally let us explain how the fornulas for the obstruction density can be generated 
without simply computing each term in Proposition ll.il This is based on the two iterative 
recursions given in |GW15 : , Lemma 2.4 and Proposition 4.9]. Firstly, given a four manifold 
with metric g and hypersurface £ labeled by a defining function s, we improve the defining 
function to 


s 



which ensures that |Vso||v = 1- It is highly propitious to work with a unit defining 
function si with |Vsi| = 1. For applications, it suffices to require |Vsi| = 1 + s £ f for 
smooth / and t of sufficiently high order. This is achieved using the recursion of [GW151 
Lemma 2.4]. For the first three examples tabulated above, si = x is already a unit 
defining function, but for the Schwarzschild case, we performed this algorithm explicitly 
using the computer software package |PM0l| . 

Supposing now that s is a unit defining function, we next improve this to a conformal 
unit defining density evaluated in the scale determined by the metric g. This can be 
achieved using the recursion of [GW151 Proposition 4.9]. In practice this simply amounts 
to repeatedly computing the leading failure of I 2 to equal unity and then adding a next 
to leading correction to a to correct for this. In Appendix [A] we present the first four 
improvement terms for a unit defining function s as well as the corresponding obstruction 
densities for hypersurface dimensions d = 1,2,3 and 4. These were computed using 
the symbolic manipulation software 4 \L VV12 . We have verified their correctness by 
checking that the resulting I 2 equals unity to the required order for explicit metric and 
hypersurface examples, including those given above. 

We note that the result for B4 quoted in the appendix, when evaluated along £ gives 
the obstruction density for five-manifold hypersurfaces expressed in terms of a unit defin¬ 
ing function. This is a novel result that is simple to use in practical applications. There 
is no in principle difficulty employing the unit defining function identities of Sections [3] 
and [4] and their higher order analogs to evaluate in terms of hypersurface invariants. 
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Appendix A. Conformal unit defining function formula and 

FIVE-MANIFOLD OBSTRUCTION DENSITY 

Let s be a defining function for a hypersurface S in a d-dimensional Riemannian 
manifold (M, g ). The following formulae improve this to Sk which determines a defining 
density cr^ in the scale g such that 

4 = 1 + a k f , k<d, 
for some smooth / € T(£M[—k]): 


s 0 = s 


si = s 0 + 

«2 = si + 

S 3 = S 2 + 


1 


2(d — 1) 

1 


s 2 V.n 


6(d — 2)(d — 1) 

1 


s 3 (2(V.n) 2 - (d - 4) V n V.n + 2 (d - 1) J) 
5d- 6 


( kW _ rj 

-j—j- (V.n) 3 - 4(2 d - 9) V.n V n V.n 


24(d-3)(d-2)(d-l) Vd- 1 
+ (d — 4)(d — 6) V 2 V.n + 3 (d — 2) A V.n + 4(2 d — 3) V.n J — 2(d — 6 )(d — 1) V^J^ 


S4 — S3 + 


1 


+ 


+ 


120(d — 4)(d — 3 )(g? — 2){d — 1) 
398d 2 — 880d— 49d 3 + 576 x2 


13 d 3 - 62 d 2 + 100d - 48 
(d — 2)(d — l) 2 


(V.n) 


(d — 2)(d — 1) 
25 d 2 - 94 d + 72 


13d 3 — 156d 2 + 524d— 384 

(V.n) 2 V n V.n H--—--V.n V 2 V.n 


d- 1 


„ lid 4 — 149d 3 + 668d 2 — 1112d + 576 ._ _ n2 

V.n AV.n +- (ri- 2)(d- !) - (V “ V ' n) 


JWI WI (3d 2 — 22d + 16)(d — 3) x2 

— (d — 4)(d — 6)(d — 8) V 3 V.n- - —--(VV.n) 2 


d- 1 


— 4(d — 4)(d — 3) AV n V.n- 3(d— 2)(d— 8) V n AV.n + 

-(V n V.n) J + 


2(13d 2 — 58d+ 72) 2 


4(4d 3 - 35d 2 + 88d - 72) , 4(d - 4)(d - 3)(d - 1) 2 


(d — 2) 

J 2 


-(V.n) 2 J 


d — 2 v y d — 2 

- 6(3d - 4)(d - 6)V.n V n J + 8(d -3)(d - 1) AJ + 2(d — 8)(d — 6)(d - 1) V 2 j) 
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When k = d, we have I 2 = 1 + cr d B^ where the obstruction density B j = Thus the 

following determine the obstruction density in dimensions d = 2,3,4 and 5. 

B\ = — (V.n ) 2 — V n V.n — J 

£?2 = — — (2AV.n + 2 V 2 V.n + 8 V.n V„V.n + 3 (V.n ) 3 + 8 V.n J + 8 V n J) 

S 3 = — TW 77 (9 V„AV.n + 12 V.n AV.n + 6 V.n V 2 V.n + 3 (VV.n ) 2 
108 V 

+ 6 (V ri V.n ) 2 + 18 (V.n ) 2 V ri V.n + 4 (V.n ) 4 
+ 9 AJ + 18 V 2 J + 36 V.n V n J + 18 V„ V.n J + 18 (V.n ) 2 j) 

B 4 = - - ( — 312V.n(V„V.n) 2 + 352V.n(V 3 V.n) — 155V.n 2 V.n 3 

46080 V 

+44 V.n 2 (V 2 V.n) — 832V.n 3 (V n V.n) + 480(V n V.n)(V 2 V.n) 

+ 96 (V 4 V.n) + 32AV 2 V.n + 288V.nAV rt V.n 
+ 256 V n AV„V.n + 96(V a V.n)V a V„V.n — 552V.n(V a V.n) V a V.n 

- 608 (V a V.?r)V„V a V.n — 1436V.n 2 AV.n — 1248(V n V.n)AV.n 

- 2176 V.nV n AV.n - 864V 2 AV.n - 288A 2 V.n - 2304 V.n(V„V.n)J 

- 832 V.n 3 J - 640(V 2 V.n)J - 384 (AV.n)J + 1024J 2 V.n + 512JV n J 

- 2496 V.n 2 V n J - 2304(V n V.n)V n J + 1536JV n J - 2432 V.nV 2 J 

- 768 V 3 J - 256AV„J - 512 (V a V.n)V a J - 2176V.nAJ - 2048V n A j) . 

The above yields B\ = B\\y, =0 and it is not difficult to show that i? 2 |s agrees with B 2 
as in Proposition 16.61 The expressions above for B 3 and B± provide an efficient way 
to compute the obstruction densities B 3 and £>4 for explicit metrics, but expressing the 
above formulae in terms of standard hypersurface invariants is rather tedious. For ( M,g ) 
flat and d = 3, this computation was performed in |GW15| and agrees with our generally 
curved expression for £>3 given in Proposition ll.il 
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